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TECHNICAL MEMORANDUM NO. 1205 


EXACT CALCULATION OF LAMINAR BOUNDARY LAYER IN 
LONGITUDINAL FLOW OVER A FLAT ELATE 
WITH HOMOGENEOUS SUCTION* 

By Rudolf Igllsch 

1. Statement of the Problem and Introduction 


Lately it has been proposed (reference 1) to reduce the friction drag 
of a body in a flow for the technically .important large Reynolds numbers 
by the following expedient: the boundary layer, normally turbulent, is 

artificially kept laminar up to high Reynolds numbers by suction. The 
reduction in friction drag thus obtained is of the order of magnitude 
of 60 to 80 percent of the turbulent friction drag, since the latter, for 
large Reynolds numbers, is several times the laminar friction drag. In 
considering the idea mentioned one has first to consider whether suction 
is a possible means of keeping the boundary layer laminar . This question 
can be answered by a theoretical investigation of the stability of the 
laminar boundary layer with suction. A knowledge, as accurate as 
possible, of the velocity distribution in the laminar boundary layer with 
suction forms the starting point for the stability investigation. 

H. Schlichting (reference 2) recently gave a survey of the present state 
of calculation of the laminar boundary layer with suction. 

The classic example for Prandtl's theory of the boundary layer 
without suction (reference 3) is the boundary layer, calculated by 
H. Blasius (reference k ) , which develops on an infinitely extended flat 
plate in longitudinal flow. Again and again one reverts for comparison 
to this example, even where the lay-out of the problems is entirely 
different. Study of the flow along the flat plate will thus have the 
same importance and significance for the boundary layer with suction or 
blowing. The developing boundary layer will, of course, no longer depend 
only on the shape of the body Immersed in the flow — a plate in the present 
study — but also on the manner of suction or blowing. The case of 

*"Exakte Berechnung der laminaren Grenzschlcht an der langs— 
angestramten ebenen Platte mit homogener Absaugung." Schriften der Deutschen 
Akademie der Luftfahrtforschung, Band 8 B, Heft 1, 1944. Report of the 
Mathematical Institute of the Technical Academy Braunschweig. Presented 
to the Germa n Academy of Avaition Research by the corresponding member 
Hermann Schlichting on January 26 , 194k. 
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continuous suction with constant suction Telocity will he of particular 
theoretical importance; it is numerically calculated in detail in the 
present report and for the general case of arbitrary suction or blowing * 

the complete system of formulas is given with which the numerical 
calculation may be performed in exactly the same manner. 

The plate is assumed to extend an infinite distance in the x direction, 
starting at the origin of coordinates (fig. 1). For negative x the flow 
is then undisturbed; thus u = TJ 0 , v = 0, where u and v signify the 

velocity components in the x and y direction, respectively, and TJ is the 
free stream velocity. The equation 


v(x,0) = v Q (x) (1) 

gives the prescribed suction law; v Q (x) < 0 signifies suction, 

v G (x) > 0 blowing. The condition u(x,0) = 0 at the wall 1 b also to 

be maintained for the boundary layer with suction or blowing. Prandtl ' s 
boundary layer equations (references 4 and 5) (with the sink effect 
neglected) read 


u ^ + v = v ' (2) 


= 0 (3) 

Sx c3y 


v signifies the kinematic viscosity. The boundary conditions (fig. 2) 


u(x,0) = 0, v(x,0) = v Q (x), u(x, ») = u 0 , u(0,y) = U 0 (4) 


also have to be satisfied. A special very simple solution of the system 
of equations (2) to (4) with the exception of the fourth boundary condition 
(4) can be given for the case of continuous homogeneous suction^ that is 
for v 0 (x) - v Q = const. < 0. In this case, according to H. Schlichting, 
(reference 6) 


V 
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u(x,y) 




v(x,y) 


(5) 


Thus, velocity distribution and hence also boundary layer thickness and 
all remaining boundary layer parameters are independent -of the distance 
along the plate x for this solution) it is cal led the asymptotic 
solution. 

For a plate provided with homogeneous suction beginning at the 
leading edge the boundary layer which starts in front with zero thickness 
must be transformed, after a certain distance, into this asymptotic 
solution. This distance in particular will be treated in detail below. 

In the first chapter (paragraphs 2 to 5) the general equations for 
the solution of the complete problem (equations (2), (3), and (L) ) are 
set upj in chapter two (paragraphs 6 to 18 ) .the case v Q (x) — v Q 

with v Q < 0, that is, the case of constant suction, is treated numerically. 

The calculation methods described here can, of course, be applied also in 
the case of a general suction law. Furthermore, the validity of these 
considerations is,, naturally, not limited to the flow over a platel; the 
numerical methods described may be used also for- airfoil profiles. 


I. THE GENERAL SYSTEM OF FORMULAS 
2 . Elimination of the Kinematic Viscosity v and 
the Free Stream Velocity U 0 


The equations (2) and (3) are simplified, if one inserts 


x = x-jVj y = JiV 


(6) 


with 


u l( z l ,7 l) = u (^y) and =v(x,y) 


■1 


This investigation has been suggested to me by Mr. H. Schlichting. 


to 
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du 

u, — — + V, 


,j_ Su-|_ 


1 Sx-L 1 ^ 3y 1 2 


^1 Q 

Sxi by ± 


according to equation (It-) tile boundary conditions 


J 


"I 


u^x^O) = 0 , v 1 (x 1 ^0) =.v q (x 1 v) =v ol (x L ), 

U^X^co) =UqJ u^O^) = U o 

also are to be satisfied. If one makes the new siibstitution 

*2 “ =7“ x l* 72 “ 7l> 

= tr- ui(xi,yi), v 2( z 2^2) = v i( x i^yi)j 
u o 

equation ( 7 ) is transformed Into 




Sx 2 ^2 dy 2 2 


^2 , 5v 2 

Sxg dy 2 


= 0 


(7) 


(8) 


(9) 


(10) 


( 11 ) 
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with the bo undar y conditions according to ( 8 ) 


( -*2 » ® ) = ® * T 2 ( ^2 ■> ) - T ol' 


UgCsgj”) = 2 , 112(0,72) “ 2 


.(“ ^ 2 )= v o 2 ( 3 C 2 )^ 


3 • Introduction of the Stream. Function ijr g as 


Independent Variable in Addition to Xg 


If one puts 


_ _ 2 . , v 

ug - - , --- 


equation (11) is aut omat ically satisfied. With 


U2(*2,72) = ^ 2 ^ 2 ^ 2 ) 


equation ( 10 ) is transformed into 


5 a . fSi + 

3xg 3^2^ \^ r 2/ 


which is equivalent to the equation 


<k_ = fa <£z_ 


(see reference 7 , p. 8 ^ 7 , formula 9 .) if one puts 

z(t 2 ,X2) = Wg 2 ^^) 


(16) 
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Because of equations (12), (13), ana (16), the following boundary conditions 
(cf. fig. 3 vhere a sink distribution is assumed) must be satisfied, aside 
from ( 15 ) : 


z(oo, 2C2) = ^ = z(Tjf 2 ,°) , z |f(x2),xg] = 0 (!7) 


with 





v o2(xq)&Z2 


(18) 


The reconversion from, the coordinates ''Kp, Xp to Xp, y P 
to (13) by 


72 



f(x2) 


U 2 (t 2 ,X2) 



dt. 


\J z(^ 2 > x 2^ 


is made according 


(19) 


(for constant Xq). 


4. Simplification of the Limits of the Eegion of Integration 


The region of integration obtains straight-line boundaries only, if 
one introduces instead of x 2 and ^ 2 tlle ne ' w ’ variables 


= t 2 - f(x2), TJ1 = F(X2) dX2 

vio 


( 20 ) 


with 


S’(x 2 ) = \|l + f , 2 (x 2 ) 


( 21 ) 


(cf. fig. 4). One may write the second equation (20) abbreviatedly in 
the form 
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^ 2 . = s(x 2 ) with the solution = S(tj^) 

(22) 

With 

z(t 2 ,x 2 ) = Z^ jT^) 

(23) 

because of 

^2- =_^- f'(xp) + — F(zo), — = — 
Sxg c&! Stii 3^2 ^1 


equation ( 15 ) becomes 



_ f * (Xg) + F(xg) = \/Z^- 

^l ^ll st x 2 

(24) 

If one designates, according to (22), 



# 

-**(:= 2 ) = “ fl [ s (t)i)] = a-(-ni) ^ 

(25) 

• 

^(x?) = f[s(ti 1 )] = 1(7^), 

( 26 ) 

equation (24) becomes 



Mu) + Mu) J- - * 

&! wii ° 6 i 

(27) 

Furthermore 

, according to ( 17 ) the boundary conditions (cf. 

fig. 5) 


Z( »^i) = 4 = Z(?!,0), Z(0,t)!) = 0 

( 28 ) 


must be satisfied. 
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5. Elimination of the Discontinuity in 
the Boundary Conditions 


We shall introduce 

6 , 

t = “r and Y = ^ 


as new variable, With. 


T(£,Y) = 


one calculates 


8Z 8 2 Z = ^1 £Z_ = _1 1 , S £ 

<n ± hi \/r 85 X 2 y 5 an! 2 ^y + &y 


(29) 


(30) 


Hence (27) becomes 




a(Y) \/Y - | b(Y)£ 


+ Yb(Y) $£■ = \JT 

8Y ^2 


(31) 


According to (28),. the boundary conditions 


T(co,Y) = k } T(0,Y) = 0 


(32) 


must be satisfied. 

Finally, a last simplification is obtainable by introduction of 


T 


= and a = \JY 


(33) 
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With 


T(r ,cr) = TCC ,Y) 


(3*0 


one calculates 


ar = Sv; _i_ Sr 
S£ St 2t j SY 


Sv l_ S^T . 1 /l S% 
Sc 2 a 5^2 hr \ T dr 2 


_1_ Sv\ 

t 2 St/ 


equation ( 31 ), vith 


&(Y)\fr = a(cr 2 )a = A(a) 


h(Y) = b(a 2 ) = B(c) 


(35) 


"becomes 


ft i f A ( 0 ) - Mid ra] + ^E(a) i S 

St 2t L 2 J 2a S a 



_i_ sA 

t 2 St/ 


or 

yfv ~ & + 2E [1 2ta ( a) + B(a)T 3 - InJt = 2T 2 crB(a) ^ 
dT 2 ^t L T - 1 Sa 

According to (32), the "boundary conditions (fig. 6 ) 

Y( “,cr) = k, V(0 ,a) = 0 


(36) 


(37) 


also must be satisfied. 

The numerical treatment of the problem (equations ( 36 ) and (37)) can 
take place in exactly the same way as that of the special case of constant 
suction velocity, treated in detail in the next chapter. 



10 


MCA 1 M Wo. 1205 


II. THE CASE OF HOMOGENEOUS SUCTION 
6. Specialization to Homogeneous Suction 


Thus we put from now on 


t 0 (i) = with v Q <; 0 
Then one obtains according to (l8) 

f(zg) = “ T o x 2 


(38) 


(39) 


according to (21) 


F(x 2 ) = 


( 40 ) 


according to (20) and (22) 


*1 = \A + V o 2 V X 2 


*1 

/l + v o 2 


(41) 


according to (25) 


a (n-, ) = v 

1 o 


(42) 


and according to (2 6) 


fc(t|i) = + \/l + v Q 2 (43) 


The differential equation ( 24 ) is therefore transformed into 


v o — + \A + v o 2 ^r~ = \!z (44) 

^l d&i 

with the boundary conditions (28) . 
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For further considerations, in order to "become independent of v Q} 

one introduces, slightly different from the procedure in the general, case 
(paragraph 5), first as new independent coordinates 


X — 


“ v o*l* 



? o 2 

ft + v 6 2 


(45) 


With 


R(X,T) = Z(S x , ni ) 


(46) 


equation (44-) then "becomes 


. (1,7) 

a » a 8 

according to (28) the "boundary conditions 

R(~,Y) = 4 = R(X,0), R(0,Y) = 0 (48) 

also must "be satisfied. 


In analogy to (29) one now uses 


£ = — and Y 

Vy 


(49) 


then according to (33) T and cr as independent variables and thus 
obtains instead of (36) th9 differential equation 


\7 


&r 2 


( 2 To + t 3 - i = 2T 2 a 
d T \ T v / Scr 


( 50 ) 


with the boundary conditions (37). 
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7. The Series Development in Powers of a 

Equation (50) is a nonlinear parabolic differential, equation of the 
second order. Usually for problems of this hind (for instance for the 
equation of heat conduction) the values, for instance, of V are 
prescribed even for 0 = 0; the reason that they are not prescribed here 
is probably that Sv/Scr contains the factor o which, for 0, tends 

toward zero. Thus 0 = 0 is a singularity of the differential equation. 
More important conclusions will result from this fact. In order to find 
out whether in spite of it the problem (equations (5°) and. (37)) can be 
regularly resolvable in the neighborhood of o — 0, one sets up the series 
development 

V(t,ct) = Y 0 (t) + Y 1 (r) a + Y 2 (t)ct 2 + . . . (51) 


There follows from it 



Comparison of coefficients of the terms not containing a in (50) gives 


\fr 0 Y 0 


1 1 


+ v 



(53) 


if the prime signifies the differentiation of Y q with respect to the 

only variable t. Aside from this nonlinear ordinary differential equation 
of the second order, according to (37) the boundary conditions 

Y o (0) = 0, V 0 («) = ^ W 


also have to be satisfied. 

For the special case v Q = 0 one has the problem of the plate in 

longitudinal flow without suction as it has been solved already by Blasius 
(reference 4) . For this problem the corresponding quantities of which 
will always be designated by an asterisk is, according to (42) and (43), 


& *(tji) = 0, *>*(tIx) = 1 


(55) 
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so that ( 27 ) Is transformed Into 

(56) 

ST)! 

If one performs instead of (lj-5) the identical transformation 

X = 6 1 , Y = tj 1 j R* = Z* (57) 

equation ( 56 ) becomes instead of (k 7 ) 

= ( 58 ) 

Furthermore one performs the same transformations as in paragraph 6 and 
obtains instead of ( 50 ) 

n/t* + 2m 43 _ 1 , 2t 2„ I H (59) 

St 2 J, V T ' 80 

with the boundary conditions (37). This problem has a solution 


y* = y*( T ) 


(60) 


which is independent of cr and determined by 


\[v* V** 1 + y**^r3 -^V^) = 0, V*(0) = 0j Y*(co) = k 


(61) 


This is the well-known solution of Blasius. 

Comparison of ( 6 l) on one hand with (53) and (5*0 on the other hand 

shows 


V q (t) = V*(t ) 


(62) 
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One knows (and can verify it Immediately by power series development 
for (6l) ) that 


\/v o ( T ) = [ = ajT + agT 2 + . . . 


(63) 


thus 


V 0 (t) = V*(t) = a^ 2 T 2 -t- 2a 1 a 2 T 3 + . . . 


(64) 


One now determines the coefficient Y 1 (t) in. the develojment (51) 
Comparison of coefficients of the power a in (50) gives 

^ v.- * i&v + ve v (* .i i^). 


or 



According to (37), because of (54), the boundary conditions 

Y 1 (0) = 0 = 7 1 («) (66) 

also have to be satisfied. If one uses for in analogy to (63) 

the expression 

\j Y l( T ) = b l T + h 2 r2 +-••! Y 1 = t i 2 ' 1 ' 2 + Sb-Lbgr 3 + . . . (67) 

comparison of coefficients for the coefficient of t necessarily gives 
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2t l 2a l + 



= 0 


since one parameter must be kept available for tie fulfilment of the second 
boundary condition (66). For the further coefficients by with v > 1 one 

obtains qualifying equations. The further coefficients V (t) with p. > 1 

in ( 51 ) also could be determined in this manner; however, this way would be 
too troublesome. We shall therefore choose another method. 

The fundamentally important result of this paragraph is the knowledge 
of V(t,0) = Y q (t) = Y*(t). According to it, any method of profile 

continuation (as it was developed for instance by L. Prandtl (reference 8), 

H. Gortler (reference 9 ) and K. Schroder (reference 10) is applicable for 
the treatment of the problem (50), (37)- However, if the expenditure of 
calculation time is not considerably increased, these methods have too 
little accuracy, particularly in the derivatives with respect to t. 

Since, in the interest of a stability investigation for not too large values 
of a, (see references 11 and 12) and, for the stability investigation of 
the present problem, (reference 13) the second derivatives of V( t, ct) 
with respect to t also are required, the application of the two methods 
to be described below is Justified: though they require more time expenditure 

than the methods described in references 8 , 9 , and 10 , they have the advantage 
of giving the derivatives of the velocity component u with respect to the 
transverse direction with a relatively high accuracy. But also aside from 
the stability investigation the present problem is of such basic signifi- 
cance that the desire for greater accuracy of the velocity distribution 
is sufficiently Justified. Inversely, the two methods in question can of 
course be used for continuation of arbitrarily prescribed velocity profiles. 


8 . A First Approximation Method 

V( T ,o) = V( t) will be approximately calculated for cbnstant a 
according to ( 50 ) under the assumption that Y( t ,cj — k) = V(t) is known; 
for a — k = 0 this assumption is correct. One replaces in (50) 
cfY/ckr approximately by the difference quotient (Y — Y)/k; if k is 
small, one-may furthermore put Y ~ Y; this will be done in ( 50 ) every- 
where Y appears in a power superior to 1. One then obtains from (50) 
the following equation for the linear approximation 

sjf Y” + (zr a + t 3 _ i. jfjyi = S^SL (y _ y) 


(68) 
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Furthermore, according to (37) the boundary conditions 

v(o) = 0, v(«) = 4 (69) 

have to he satisfied. The considerations of paragraph 7 lead to the series 
development in powers 

]jv( t) = c x T + c 2 r 2 + C^t3 4 . + . . . ( 70 ) 

thus 

T(T) = dgjT 2 + d 3 r3 + d^ + d^T 5 + . . . ( 7 1) 


assuming as analogous 

= a^r + • « « + aj|T^" + . 

T = hgT 2 + . . . + h^T^ + 



(72) 


in the case a — k = 0 this expression according to (63) and (64) is 
Justified. The ay and b v are connected by the formulas 


*4 

b 5 


a 2 
a i 

a i 

2 a l a 2 


2 a l a 3 + 

a 3 


a^ 


/*2 
b 3 :2a ± 

(b^ - ag 2 ) :2& 1 
(b^ — 2aga 3 ) :2a^ 


> 


( 73 ) 
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* 


Comparison with the coefficients in ( 68 ) yields for calculation of 
the d v the formulas 


T 2 


2 x 3a^d^ + 1 x 2 agd 2 + 2cL^^2a — a^) — 3d^a^ = 0 


t3 


3 X 4aj_dj 4 . +■ 2 X 3®^d^ + 1 X 2a^dg — 2a^dg + 3 ^3 ( ^ cr — ag) 


— Ua^djj. = 0 



k X 5aid^ + 3 X ta^d^ + 2 X Ba^d^ + 1 X 2a^dg + 2d2(l — a^) 


— Ba^d^ + lkLi,.(2<x — a^) — 5apd^ = (dg — b 2 )cr/k 


Further calculation yields the following formulas for. determination 
of ^ : 


d 3 


ka&2 

3aj_ 


dli- =“ 


Bd^Cag + 2a) 


8&i 


d 5 - 


Sd^Cag + a) - 3d 3 a 3 - 2^(1 - g/k) - 2bgg/k ^ 

15&! 


y 
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The further numerical treatment will he discussed in paragraph 12. 


9. Yariables Suitable for Tabulation 

We shall now consider what dimensionless Yariables are most suitable 
as characteristics and what connection they hhYe with the coordinates T 
and cr forming the basis of our calculation. We shall in 


V(t ,a) = 


( 75 ) 


again introduce the original Yariables. According to (33), (49), and (45) 


thus 


N( T,cr) 


u(i/t,\|y) = u 







(76) 


Further one obtains according to (41), ( 39 ), anil (20) 


U(T, ff ) 


U 


2 + V o x 2 

V & 



( 77 ) 


From. (19) one concludes with (39) 
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If one, furthermore, takes (6) into consideration, one has between the 
original coordinates x,y and the end coordinates T jC r the connection 
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are therefore advisable as dimensionless variables. (Cq = — v 0 /U 0 ) is the 
so-called mass coefficient of the suction, defined by Q = —v 0 bZ = CQU 0 bZ, 

with Q signifying the suction quantity, b and l vidth and length, 
respectively, of the plate section considered. 
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tlms with (85) 


— v_S* 




— -^Ady* = a 




the momentum thickness 


— /i - — W 

u 0 v u 0 ; 


J poo p CO 

[tf-rj— lE-^l 


the form parameter 


5* _ "^oS* ~V 


and, finally, the wall shearing stress (p = density, pv = h) 


-■© 


dy/yi=0 


T o 8 * -v Q 8 * 


lie; 


ly *=0 
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11. The Blasius Solution without Suction 

Since the Blasius solution forms the starting point, for the .entire 
calculation, we shall first sketch its determination in our coor dinat es. 
For the numerical calculation of the solution V*(t) of the problem ( 6 l) 
the method of Adams was used, and in particular the extrapolation method 
(cf. G. Schulz (reference 14), Nr. 108) . For practical calculation with 
the calculating machine, the following formula (symbols as in Schulz* 
report) is more convenient than the formula ( 2 ) (cf. the equation (l) 
on p.. 116, Schulz) which uses the differences of the function f(x,y) 


sv+i - y v ° Ms. 291667$, - 2.458333$, 

+ 1. 541667$, _£ - 0 . 375$, _ 3 ) 


(93) 


equation ( 93 ) operates only with the values of the function f(x,y) instead 
of using their differences. This formula was applied with the width of 
interval h = 0.05 to the system of differential equations following 
from ( 6 l ) 


W** 



= W* 


(9 1 0 


The necessary four initial values were determined by series development. 

If one uses for \Jv* and V* the expression ( 72 ), the insult obtained 
by comparison of coefficients for ( 6 l) is 


a 2 ' a 3 - a 5 - a 6 ’ °> \ ■ T5 

bg = M = ^5 = = 0 , \>2 = — 2 a^a^ 


(95) 


For determination of the as yet undetermined coefficient a^ or \>2 the 
value 


f D * *( 0 ) = 1.328242 


(96) 
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was used wiiicli can "be found in reference 15, P. 551* Accord ing to 
Hcwarthj p. 550, formula (2), the variable is 


in our symbols 



(97) 


With v Q * 0 equation (76) is, when (57) is taken into aonsideration, 
transformed into 


U*(T) 



(98) 


equation (82) becomes 


•Ul*(xi,yi) = 5a U*(t) 


(99) 


and (84) 


= 7 - v 7i = Vy 2 = 2V 


n t 


nr 


-iff- = 2 . [EE 

**(r) \u o u 


TdT 

U*(T) 


( 100 ) 


Hence (97) is transformed into 


PT 




0 


TdT 

tj*(t) 


( 101 ) 


According to Howarth (reference 3.5, p. 548) f’^ij) is identical with our 
velocity U*^ though as function of q. Thus one has 


...(,) .SEClE 


dT dtj 
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or, talcing (101) into consideration. 




dT 


dT 


U*(t) 


(102) 


If one puts here T = tj 
consideration 


= 0, one obtains, talcing the development (72) into 


ap 2 = -b 2 = V r (°W2* 


(103) 


For improvement of the values attained by means of Adams’ extrapolation 
formula (93) the interpolation formula of Adam (Schulz, reference 14, 
p. 121, last formula (8)) was used which, when the differences are eliminated, 
appears as follows: 


y. 


(k+1) 


m 


- y. 


m-l 


(k+1) _ 


= h 


0 . 375 f m ^ + 0.791667^3^ 


- 0. 208334^2 + O.o 4 l 667f m _ 3 ( k ) 


3 


(104) 


The numerical values thus calculated of 

-H-, i. and U 2ElL 

"o a,* H 0 V 2 


are compiled at the beginning of table 1 as function of 


T)* 



The characteristic boundary layer parameters are to be defined as 
in 10. If one introduces as dimensionless variable (cf. also (100)) 
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r\* = y, 


Ha .^[2 


n t 


VX 


TdT 
U*( T) 
2 


( 105 ) 


(this is the quantity 2 r\ used by Howarth, cf. (97) and (101)), there 
follows according to (86) 


5 * 



/**> eo 



O oo 



(106) 


In analogous manner one obtains, according to (88), for the momentum 
thickness 


r>°° 



( 107 ) 


The form parameter again is 



and the wall shearing stress results, according to (91), as 


l£^ = 5* 




(108) 


( 109 ) 


The following numerical values result; 
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• 

&*\l 

fu - 

= 1.7207, 

— ££ = 0.571k 


V 

v x 

n u 0 



Ha = 0.66U1, 

= 2.591 


\ 

VX 

« J 


( 110 ) 


12. The Numerical Determination of the Approximate 

Solution V( Tjo) According to the First Method 

First, a solution of (68) and (69) was determined for o = 0.1 hence, 
one had to put 


T( T) = v(r ,0) = V*( t) 


The series development described in paragraph 8 can "be performed according 
to the formulas given there, only the coefficient dg (and hence c-j.) 

remains undetermined. Thus one has to assume first an arbitrary value 
of d2 and has to determine it in the end from the second condition (69). 

It will be usefull to choose d^ ^ bg . The calculation was made in a manner 

similar to one applied by Hans Joachim Luckert (reference 16). Two 
different values b g (let us assume b 21 and b 22 ) were selected, the 

corresponding values of V( t) for r =0; 0.05; 0.1; 0.15 were 
calculated according to the series (71)* then the two solutions, to be 
called T^( T ) and Vg(T)^ were followed up further by means of Adams' 

method as in paragraph 11 (cf. (93) and (10k)); the two values 
and b 22 were selected in such a manner that V^(t) > V(t) 
and V 2 (t) < Y(t) are valid for small values of t. As soon 
as V 1 ( T ) or V 2 ( T ) a Point T 0 succeeded the function V(t) or 
deviated from it by more than 1, two new functions V^t) and V^(t) 
were introduced by 


V 3 (t) = mT^T) + nV 2 (T), V 4 (t) = hY^t) + T 2 (t) 


( 111 ) 
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with 


r. W-W jn . l-M 
Vo> - Vo> 

Vo) - Vo) 


( 112 ) 


( T 0 ) > Y(t q ) and T^(t q ) < V(t q ) were selected so that behind t q 
again V 3 ( T ) lies above, Y^(t) below Y(t). After calculation of the. 
starting values of Y^ and required for the Adams method and of the 

corresponding derivatives at the points t q j — 0.05; t q — 0.1; 
and t q - 0.15 according to (ill), Y^ and maybe calculated 

further according to Adams. Of course Y^( T 0 ) = Y 1 ( t q ) or Y^(t q ) = V p (r q ) 
may be selected if that seems suitable. Possibly the pair Y^ ( t ) , Y^(t) 
also must be replaced in the same manner by Y,_(t), Y^(t), etc. The 
calculation was carried through up to a pair of functions Y pr[ _ n ( t ) , Y ?n (T) 
until it became evident that for instance Y^^t) with the assumption 
Ypn_i ( T i ) = 4 -did not any longer show a tendency .to deviate from the 
value 4 to a noteworthy degree. This was shown to occur at = 3; 

for the larger 0 values this point shifted to the left (cf. paragraph l8) . 
The desired solution of (68) for 0 = 0.1 then is 

Y(Tj 0.1; 0.1) = Yg^CT) (113) 


This solution was calculated for all t values, taking the formulas (111) 
and (112) into consideration. 

However, the accuracy of this Y( T ; 0.1; 0.1) is not yet sufficient. 
Thus Y(t ; 0.1; 0.05) also was calculated, by subdividing the 
Interval 0=0 to 0 = 0.1 into two partial intervals of half the 
width; the numerical treatment was the same. 

In order to obtain an estimate of the accuracy, the 0 interval of 
the width 0.1 finally was divided into four equal parts and 
thus V( T ; 0.1; 0.025) obtained. By linear extrapolation one obtained 
from 
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Y(t ; 0.1; 0.1) and Y(t; 0.1; 0.05) 


Y^( T ; 0.1) = 2Y( T ; 0.1; 0.05) -Y<T; 0.1; 0.1) (lU) 


furthermore by quadratic extrapolation from. 


Y(T; 0.1; 0.025 ) s V(T ; 0.1; 0.05) and Y(T; 0.1; 0.1) 


V (t- 0.1) = 2.66666TY(t ; o.l; 0.025) 


(115) 


-2Y( T ; 0.1; 0 . 05 ) + 0.333333V( T J 0.1; 0.1) 


The functions (1120 and (115) differed at most by one unit of the third 
decimal, their first and second derivatives with respect to T at most 
by 2 or 6 units of the third decimal. Finally one put 


Y(T; 0.1) = Y 2 (T; 0.1) 


(116) 


In the same way one calculated 


Y( t ,<t) for -a = 0.2; 0.3; . . . ; 1.0 


using 


V(T) = Y(T; a - 0.1) 


etc.; however, one always puts 


Y(t,cf) = Yq( T; a) Tor cr > 0.1 


( 11 ?) 
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The sub-division of the a interval of the width 0.1 into four equal 
parts was not undertaken. Actually, (verified hy the sketches of the 
curves), the variation of V(?,.c) and its derivatives as well as the 
velocity variation calculated from it with its derivatives becomes 
more and more uniform with increasing a values. 


13. A Second Approximation Method 

Following, another approximation method will be described which 
yields results much more rapidly; it requires only about one third of 
the time expenditure of the former method and operates with at least 
the same accuracy. The problem to be solved is (50), (37). By 
interpolation one dete rmine s 


— i [-0.33333 V(t , cr - 3k) + 1.5Y(r,c - 2k) 
da k L 

(117) 

- 3V(t ,a - k) + 1.83333 V(t ,a)J 


V( +,a) is calculated under the assumption that all V(T,a — nk) for 
n = 1, 2, 3, k are known. (Actually this method could, therefore, have 
been used already in the calculation of V(Tj OA), For V(r,a) one 
first determined by extrapolation the approximate value 


V*( T ,a) =-V(T,cr _ lvk) + 4 V(t , a - 3k) - 6v(r,a - 2k) 


+ 4-Y(t,ct - k) 


(118) 


The (50) is approximately replaced by 


d% + 

— + 


d 7 / 2ra + t 3 _ A = 2 T g a cY 
dr \ \[ y * V da 


(119) 


If one puts according to (117) 
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g (T) = 0.33333 T(t ,a - 3k) + 1.5 V(t,ct - 2k) + 3 V(t ,o - k) (120) 

i/V*k L J 


equation (119 ) is transformed, into 


& + |t/ 21^ _ _ ,.83333 2 ^ y , g ( T , 

^ a T l Jv* -rj k\/Y* 


( 121 ) 


This equation is solved in the manner described in paragraph 12 with the 
boundary conditions 


V(0) = 0; V(oo) = 4 


( 122 ) 


Split into a system of differential equations of the first order, it will 
appear as follows: 


W 


w* =-w 


2T q J_ r3 -!)+ 1.83333 — | V + g(T) 
tfV* T / kVv* 


(123) 


One also could have determined 
to ( 117 ) by interpolation from 


dV( T ,cr)/ck7 instead of according 


V(T,q - k), T(t , a - 2k), T(t , a - 3k), T(r,a - 4k) 


by extrapolation. In place of (121) one would then have obtained 
for V( T ,cr) a linear differential equation of the second order in which 
the term with V is missing and which thus can be solved by quadratures 
( Simpson* s rule). This method may perhaps seem more convenient to many 
readers 
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For determination of the starting values for Adams’ method one again 
sets up the series develojments (70), (71) and analogously 

][f* = c 1 *r*+ c 2 *t^ + . . . + C5* t 5 + . . . (124) 

v* = a^T 2 + . . . + d 6 *T 6 + . . . (125) 

According to (ll8) 

d v * = -dv(cr — 4k) + 4dy(cr — 3k) — 6(^(0 — 2k) + 4d v (cf — k) (12o) 
Furthermore according to (120) 

\/T*g(T) = + . . . (127) 

with 

A 4 - 20 eE 0.33333* 2 (ct - 3k) + 1.5b 2 (cr - 2k) - 3b 2 <cT - k)| (128) 

The comparison of coefficients is made, instead of with (121), to tetter 
advantage with 

\fr* ^ (2 Ter + T 3 - i \Jv*\ - 1.83333 T =#*g( T ) (129) 

8+^ S T \ J k 

■One obtains 
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4 3 (<t) = - 


4gdg(g) 

3c x * 


<^(g) 


( 3°g* + 6g)d-3(g) 


8 ct* 


^(g) = 


Ajj. -{?. - 3*66667 j^td^g) + 303*43(0) + 8(g + cg*)di,.(g) 


(130) 




150 ^ 




It is useful to select for the arbitrary (^(s) 


dgCa)** dg* 


(131) 


The first calculation according to this method was carried out for g = 1.2 
and k = 0.2; the maximum difference between Y and Y* amounted to less 
than 5 units of the fourth decimal. Nevertheless one must not think the 
application of Adams’ method superfluous; it should be recalled that one 
obtains in this manner 


SY /5 t and S 2 Y/&t? 


with appropriate accuracy which is not well possible by means of extrapolation. 
And some time, of course, this good agreement between V and Y* does 
have to stop. For g = 1.6, for instance, the difference was already one 
unit of the third decimal. 

The following g values were treated according to the method 
described above: 


g = 1.2; 1.4; 1.6; 2.0; 2.4; 3«2; 4.0 
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lit-. The Calculation of the Reduced. Telocity 
and Its Derivatives 

According to (75) ( 82 ) one finds immediately 


£--i i/w = 


U- 


(132) 


first as a function of T . The quantity y* desired as independent 
variable instead of T was calculated according to ( 85 ) by means of 
Simpson's rule. Furthermore one obtains, again by means of ( 85 ), 


JL = A. 3a £L_ = 1 VfT ) _ \JyCjI _ w&l 

XJ 0 By* n o St By* 2 2\/V(t) 2ot 8gt 


and according to ( 71 ) 


1 Bu _ 

U By* 4c 
o 


for t = 0 


(133) 


(134) 


In the same way one calculates 

_1_ = TV«(r) — W(t) Mr) = , 1 (TVf' -W)i/7 

U o By* 2 8 oT 2 2cjt 16 c 2 t 3 


and 


1 B 2 u 
U o By* 2 



for T = 0 


(135) 


(136) 


The quantity V (t ) appearing in (135) was already determined in the 
calculation according to Adams' method. 

Table 1 added, for the various values of 


I = 1/2 a 2 



34 


NACA IM No. 1205 


the calculated numerical values of 


u/u 0 i u o 1 du/ 5 y* 


and U 0 u/Sy* 2 


as a function of y* including the Blasius profile for | = 0 and the 
asymptotic suction profile (cf. (5)). Concerning the accuracy of all 
tables it should "be noted that the last decimal given cannot always he 
fully guaranteed; particularly for smaller values of the independent 
variable small inaccuracies in the last decimal may occur^. In the 
following way one obtains another estimate of the accuracy of the 
numerical material. Fulfilment of Prandtl*s boundary layer equation ( 2 ) 
for y = 0, that is at the plate , yields the result that there muBt be 


1 5 ^u 1_ chi 

*o dy* 2 U o Sy* 


7 * 


0 


In the approximate calculation this condition was not used. 
Fig. 7 shows a sketch of u/U 0 as a function of 


y* = 


The velocity distributions approach for large | the asymptotic suction 
profile (£=")• In fig. 8 the curves u/U Q against y/&* are 

illustrated with a shifted starting point. One can see that these curves 
hardly differ in shape. The curves of 


S*/U 0 1 du/dy 


illustrated In fig. 9 In the same manner, have a slightly more varied 
aspect, whereas 


^u/dy 2 


^The numerical calculations were verified by Mr. H. Schaefer on the 
basis of the so-called momentum control 


r o(*) 


= U 


2 dS 


dx V ° 


JL 


It proved useful to integrate this equation in the x direction. Agreement 
according to the accuracy Indicated above resulted. 
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(fig. 10) for- the different values of the parameter | shows an entirely 
different shape. 

The characteristic boundary layer parameters 


— v 0 &*/v and — v 0 t 3 /v 


result from (87) and (89) , respectively, by means of Simpson's rule. 
Thereafter one may calculate immediately the form parameter 6*/8 
according, to (90) and, according to (92), easily 


t 0 s*/hu 0 


With (85) taken into consideration, there follows immediately from (87) 
and (106) or (89) and (107), respectively. 



( 137 ) 


and 



d\flF 


a=0 



o.66ia 


(138) 


the numerical values ■were taken from (110). 

The four characteristic constants are to be found in table 2 as 
function of I and VT , respectively. Fig. 11 gives an illustration 
by means of a sketch; the tangents at the starting point according to 
equations (137) and- (138) are plotted there. 


15. The Stream Line Pattern 

As stream function we shall introduce the stream function 1 Hl,y*) 
with the use of the dimensionless coordinates; it is defined by 
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St _ v _ u_ 

~ y 0 ’ by* U Q 


(139) 


For that purpose one puts 


♦ - | (-^ o)*2 


(UO) 


Actually, according to (83) and (85), 


£ = \ i- y o)^2 


(ito) 


hence according to (llt-O), (13), and (9), 

^ = i (_ v ) ^2 2 y 

Si 2 bxQ (_ Tq ) 2 -v o 

furthermore according to (8A) and (85) 


y* = (“ y o) y 2 


(1^2) 


thus according to (140), (13), and (8) 


ZL . 1 ( _ v ) _j=_ = u. ■ 

Sy* 2 0 Sy 2 - v D TJ 0 

The numerical calculation of the function + may he carried out as 
follows: From (79) one concludes 


*2 = V^t 2 -^0*2 


with (83) there results 
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thus according to (lhO) 


ijr = 1 ct(t 2 + a) 


(1^3) 


This value corresponds to 


| = i a 2 and y*(T,c) 


(1M) 


the latter has to he calculated according to (85). 

According to the tables calculated in this way the pattern of the 
stream lines ^ = const, was constructed by means of linear interpolation. 
It must be noted that formula (lV3) yields for | = 0 and hence a = 0 ^ 
only the singular stream line (one point) ijr = 0 with y* = 0 (cf. (8->)). 
Since , however, for 5 = 0 , u/U 0 must equal 1 , one has because of the 
second equation (139) 


(0,y*) = J* 


(lh5) 


This result is naturally yielded also by our system of formulas: 
according to (85), namely. 


P T 




a 



TdT 

U(T,c) 


(l b6) 


whereas from (1^3) follows 



( 1 ^ 7 ) 


The right side of (lh6) goes, due to the monotone growth of U( t,c) as 
function of t from 0 to 2, with increasing t from higher values 
toward 
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2 


TT 

A) 



T 2 


2 


so that for 0-^0 and large t (and only these have to he considered 
for 0 -» 0 , according to ( 1 * 1 - 3 )*} the right sides of (l46) and (l*t 7 ) agree. 
That, however, leads to ( 1 * 1 - 5 ). 

Fig. 12 presents the stream line pattern, fig. 13 a partial section 
for I— values from 0 to 1 in magnified scale. Particularly noteworthy 
is the steepening ‘of the stream lines at the beginning of the plate 1=0 
where all stream lines have a vertical tangent (cf. also paragraph 1 6 ). 

All stream lines approach the wall vertically, since there v = v Q 

and u = 0 . 


1 6 . The Velocity Component Perpendicular to the Wall 
The velocity component 


v (£ *y*)/-v Q 


perpendicular to the wall results from ( 139 ). On© calculates 


Sty _ 3^ 3cr 3jf 3 t 

as " 30 ai + st si 


Now, according to ( 85 ), 


dl = 0 d 0 , dy* 


20 


U(r,0) 


, v <V* , 

d*r + < 3.0 

30 


(148) 


( 149 ) 


with 


3y* _ y* 
30 0 


- 20 


r\ 


Sglli£l 4T 

3c 


Jo 


U 2 (t,0) 


(150) 


By further calculation one obtains from (1*1-9) 
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dcM., dT , SLV1 /W* _ SC 

2a T 8c a/ 


Since from (1^3) follows 


S£ 

Sc 


1 (** ♦ 2.). * 

2 i T 


= cyT 


one obtains, according to (1^8), 


_ v_ = 2S = i. ( t 2 + 2s ) _ StllSl M 
v o 8^ 2c 2c Sc 


(151) 


(152) 


(153) 


This is the formula for calculation of tbe velocity component perpendicular 
to tbe wall. Sy*/Sc can be determined by interpolation, Bince y*( T ,c) 
was required anyway for tbe former calculations. 

For T = 0, where 


U(0,c) = 0 and Sy*(0,c)/Sc = 0 


necessarily follows from (153) r /~ r Q = —I 1, that is v = v q . For r ^ 0 
and a -> 0 tbe decisive factor is, according to ( 153 ) and- (150), 


llm 

c 0 a \ 2 2 a / 


lim - \- t 2 
cr 0 y 2 



(15^) 


equation (85) was taken into consideration. Due to tbe monotone growth 
of U( T ,e) from 0 to 2 with t the bracket is negative so that one 
obtains 


v 


— v 


o 


+ for t =0 and c = 0 


(155) 


This, however, means nothing else but the steepening of the stream lines 
for I =0 pointed out at the end of the preceding paragraph. Moreover, 
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one may conclude from (154) in connection with the two equations ( 85 ), 
that for g 0 and y* different- from zero 


— = — lim U(t,o) — — -t 2 + . . . 

“ T o \/ 2 | t ® - U(T , 0 ) 2 

LUO J 


Instead of calculating lim T “ it is sufficient to substitute that 
value t* for which TJ( t*, 0 ) =2 is satisfied with adequate accuracy. 
According to the calculations described in paragraph 11 one can 
put t* = 3 . Then one obtains 


- T o ^ 


U(t,0) 


- 4.5 + 


l/§f 4.Jj 


y* being the value defined by (105). The numerical calculation gives 
B7 

the result 


v 0.8448 




( 156 ) 


For performing the calculation of the 8 y*/ 8 c for t > 3 appearing 
in ( 153 ) one should further note that, according to ( 85 ), 


y*( T ) = y*(3) + | ( t2 - 9) (t ^ 3) 


(157) 


Hence follows immediately 


8y*(T) _ 8y*(3) + 1 #2 

8 c 8 c 2 


-9) (t \ 3) 


( 158 ) 
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hi 


and thus from. (153) that 


v(t) _ v(3) 


for t > 3 




— V, 


(159) 


In fig. lh the values of 
the various values of 


v/v Q are plotted as function of y* for 


17 . The Drag of the Plate 


The total friction drag for the plate of the width b 
length l wetted on one side is 


and the 


11 

W = b T 0 (x)dx 
00 


(160) 


If one defines in the customary manner the coefficient of the friction 
drag by 



(161) 


and puts 

T q8* 

f( |) = Io_ m , ^ U P „ T ° 

T Oo3 v O$* p(— v„) Uo 

V 


(162) 
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(regarding the value T Op0 for the "asymptotic solution" (reference 6) 
(also see paragraph 18) c f 160013188 "with 


«1 


V \ 2 UpZ c gU Q Z 

v v Q v 


(163) 


ns 


c f =2^-i 

u 0 h 


Uo 


Is- dl = 2 -Is. 

o CO 


Uo 


F(Si) = 2 cqF(! I ) (164) 


with 


ps 


F(4 2 ) = 


*7 


f(4 )d| 


(165) 


UO 


One concludes immediately; 

For 4 , -> “ one obtains 
i ♦ 


F(4 Z ) 1, 



C f oo— 2Cq 


for t z 0 one obtains 


F(4 ) 4 Q-iffljj:; c ->^0 1,3282 . 1.3282 




/h 



The first is the drag law of the asymptotic solution, the latter Blasius’ 
drag law for the case without suction. The drag law in the entrance 
region is given by equations (163) and (164). The drag coefficients for 
various mass coefficients Cq coincide in one curve whon 
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c f /c fM = c f / 2 c Q = F(l 2 ) 


* 

is plotted against 


c q % 0 2A = £ z 


(Fig. 15). 

For calculation of the integral (l 6 o) one introduces according 
to ( 85 ) a as integration variable : 






Jo 


f^ajadc 


with 


2 a l~ = *1 and f l(’ a ^ = f (^) 


(166) 


(167) 


According to ( 87 ) and (92) as well as to (106), (109) and (110) 
5*v n \ T _ 5 * 

= 1.7207a . . ., — — = 0.5714 + . . . 

V V2 

t 

hence 

*(!,) .-S- H v5M2±da + . . .■.i 26 2 1 2Et + . . . 
a 2 2 Jo 1.7207 ^2 1.7207 


= 0.6641 -4=. 

✓*T 


(168) 
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in agreement with the statement made in connection with ( 165 ). The integral 
appearing in (l 66 ) is a definite integral and was evaluated according to 
Simpson’s rule; first with the width of interval h = 0.2; thus one 
obtained the values of Ffl^) for = 0, 0.2, 0.4, . . . From them the 

value of F for 0 ^ = 0.1 was determined according to Newton's Inter^- 
polation formula; by repeated successive application of Simpson's rule 
the F values for = 0.3, 0.5, 0.7, . . . were obtained. 

The last column of table 2 gives the values of the function F(t,). 

3 * 

Fig. 15 shows a sketch of this function. In fig. 1 6 the curves 10 J c^ 
are plotted as function of TJ 0 z/v with cq = — v 0 /U 0 as parameter; for 

comparison the curves are drawn as dashed lines which illustrate for the 
flow without suction the transition from laminar to* turbulent flow or, 
respectively, the fully turbulent flow (reference 1 ?). 

The difference between the fully turbulent drag curve and a 
curve cq = const, gives the drag reduction by means of keeping the 

boundary layer laminar, under the presupposition that the laminar boundary 
layer is stable for the respective mass coefficient, that is, that no 
transition to the turbulent flow type occurs^. 


18 . The Asymptotic Solution 

According to Schlichting (reference 6 ) the problem ( 2 ), ( 3 ) is solved 
by the expression 


v = v Q , u = U Q (1 - e-T*) 


(169) 


lt lnformation on the mass coefficients required for maintenance of a 
laminar boundary layer Is given in an Investigation by A. Ulrich 

(see reference 13 ) in which the laminar velocity profiles calculated 
in the present report were investigated for stability. The mass 


coefficient 


= 1.2 X 10 


,—4 


was found to be sufficient for the 


Qcrit 

maintenance of a laminar boundary layer. Hence reductions in friction 
drag of the order of magnitude of 'JO to 80 percent of the fully turbulent 
friction drag result for the Reynolds number in the region 


U l ft = 5 X 10° to 10 


8 


Important In practice. 
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■which also satisfies the boundary conditions (^) with the exception of 
the fourth, the latter only for y*->e°. With increasing y* the 
individual velocity profiles u/U Q ‘will tend toward this "asymptotic 

solution" ( 169 ). We shall reconvert ( 169 ) to our calculation 
coordinates 0 and T : 



One takes ( 85 ) and (132) into consideration: 

nr 



Differentiation with respect to T gives 


SaL.liidL 

Vv 1 _ 1 


b ft t-ft 

= — 20TdT 

1 -iv? 

2 


and integration 


\/V + 2Zn ( 1 - i- 

V 2 




(170) 


(171) 


The additive integration constant which actually should appear on the 
right must disappear since according to ( 170 ) for r = 0 v' v also must 
disappear . 

From (171) one concludes for instance: For constant r 'S V tends 
with growing 0 increasingly toward 2 j this is the reason for the remark 
tfiade in paragraph 12 that the quantity called t there decreases with 
growing 0. 

Equation (171) assumes a clearer form if one introduces according 
to (33), ( 3*0 and ( 4-9) Z and Y as independent variables: 



46 


MCA TM Ho. 1205 


\/H + 2ln f 1 - i Vf J = -Z (172) 


This could, hare been confirmed, directly from (47). 

In all figures the values referred to the' asymptotic solution which 
can easily he calculated are marked, in. It is shown that for a = 4 all 
flow characteristics lie sufficiently close to the corresponding values of 
the asymptotic solution. 


SUMMARY 


A general calculation method is given with which to determine the 
boundary layer developing along an infinitely extended flat plate under 
the influence of an arbitrary suction or blowing law. For the special 
case of homogeneous suction the numerical calculation of the velocity 
within the boundary layer is performed completely. Two numerical methods 
are applied which are both fundamentally based on the approximate solution 
of an ordinary linear differential equation of the second order for the 
calculation of each velocity profile vertical to the plate; thus one 
attains a sufficiently accurate determination of the first and second 
derivatives of the velocity vertical to the wall as well. Both methods . 
can also be applied for the continuation of an arbitrarily prescribed 
velocity profile. * 

Translated by Mary L. Mahler, 

National Advisory Committee 
for Aeronautics 
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TABLE I 


THE TELOCITY DISTRIBUTION WITH ITS FIRST AND 


SECOND DERIVATIVE FOR VARIOUS I 


B 

sss 

= 0 Blaslus 

cr 

II 

0 

H 

IAt 

= 0.005 


■ 

u 

1 8u 

1 ci 2 u 


u 


1 c) 2 !! 

■1 

Uo 

u 0 dn* 


y* - — v" 

Xio 



0 

0 

0.332 

0 

0 

0 

5.322 

-5.329 

0 .2064 

0 .0685 

0.332 

-0 .001 

0 .0137 

0 .0726 

5.248 

-5 .481 

C .4128 

0 .1370 

0.331 

-0 .005 

0 .0275 

0 .1445 

5.169 

-6 .100 

0 .6193 

0 .2053 

0.330 

-0 .011 

0 .0414 

0 .2156 

5.077 

-7 .149 

0-7935 

0.2733 

0.327 

-0 .018 

0 .0554 

0.2858 

4.968 

-8.595 

1 .0007 

0.3406 

0.322 

-0 .028 

0 .0694 

0.3546 

4.835 

-10.389 

1.2088 

0 .4069 

0.315 

-0 .o4o 

0 .0836 

0.4220 

4.673 

-12 .459 

1 .4179 

0 .4718 

0.305 

-0 .053 

0 .0979 

0 .4874 

4.479 

-14 .705 

1 .6286 

0 .5348 

0.293 

-0 .066 

0.1123 

0 .5505 

4.250 

-17 .002 

1.8412 

0.5955 

0.277 

-0 .079 

0 .1269 

0 .6108 

3.985 

-19.204 

2 .0563 

0 .6532 

0.259 

-0 .092 

0 .1418 

0 .6678 

3.685 

-21 .150 

2 .2744 

0 .7074 

0.238 

-0 .102 

0 .1569 

0 .7210 

3.353 

-22 .676 

2.4964 

0 .7575 

0.214 

-0 .109 

0 .1723 

0.7699 

2.995 

-23 .626 

2 .7227 

0.8032 

0.189 

-0 .113 

0 .1881 

0.8142 

2.620 

-23.880' 

2.9544 


0.163 

-0 .113 

0.2043 

0.8535 

2.237 

-23.361 

3 .1922 

1 

0.136 

-0 .109 

0 .2209 

0.8875 

1.858 

-22 .062 

3 .4371 

0.9097 

0.111 

-0 .100 

0 .2381 

0 .9163 

1.495 

-20 .051 

3.69OO 

0 .9347 

0 .087 

-0 .089 

0 .2558 

0.9398 

1.161 

-17 -472 

3.9519 

0 .9545 

0.065 

-0 .075 

0 .2743 

0.9584 

0 .866 

-14-537 

4 .2237 

0 .9697 

0.047 

-0 .060 

0 .2934 

0.9725 

0.617 

-11 .491 

4.5063 

0 .9808 

0.032 

-0 .045 

0.3134 

0.9827 

0.425 

-8.586 

4.8007 

0.9885 

0.021 

-0 .032 

0.3341 

0 .9898 

0 .267 

-6 .027 

5.1074 

0 .9935 

0.013 

-0 .022 

0.3558 

O'. 9943 

0 .160 

-3.949 

5 .4271 

0 .9966 

0.007 

-0 .013 

; 0.3784 

0.9970 

0 .089 

-2.399 

5.7603 

0 .9984 

o.oo4 

-0 .008 

i 0 .4020 

0 .9986 

0.046 

-1.342 

6 .1072 

0 .9993 

0 .002 

-0 .oo4 

j 0 .4265 

0 .9994 

0.022 

-0.690 

6 .468o 

0 .9997 

0.001 

-0 .002 

0.4520 

0.9998 

■ 0 .009 

-0.323 

6.8428 

0 .9999 

0 

-0 .001 

0.4785 

0.9999 

o.oo4 

-0 .136 

7.2317 

1.0000 

0 

0 

0 .5060 

1 

0.001 

-0 .051 

7-6348 

1 .0000 

0 

0 

0 .5344 

1 

0 

-0.012 

8.0519 

1 .0000 

0 

0 

0 .5640 

1 

0 

-0 .006 
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TEE "VELOCITY DISTRIBUTION "WITH ITS EIRST AND 


SECOND DERIVATIVE FOR VARIOUS £ 



% = 0.02 


1 Su 1 5 2 u 
Uq ciy* U 0 by*£ 


0 

0 .0260 
0 .0522 
O.O787 

0 .1054 
0 .1324 
0.1597 
0 .1871+ 
0 .2151+ 
0 .2439 

0.2730 
0.3026 
0.3329 
0.3640 
0.3960 
0.4290 
0 .4631 
0.4984 
0.5351 
0.5733 
0 .6131 
0 .6546 

0 .6979 
0 .7431 

0 .7902 

0.8392 
0.8902 
0 .9432 
0.9982 

1 o0552 
1 .1142 


0 

0 .0766 

0 .1519 
0.2257 
0.2980 
0 .3684 

0 .4310 

0 .5027 
0 .5658 
0 .6257 
0.6820 
0 .7342 
0 .7820 
0.8249 
0 .8627 
0 .8952 
0 .9225 

0 .9442 
0 .9621 
0 .9751 

0.9845 

0 .9909 

0 .9950 

0.9974 
0 .9988 

0 .9995 
0 .9998 
0 .9999 

1 
1 
1 


2.986 

2.910 

2.831 

2.748 

2.658 

2.558 

2.445 

2.319 

2.178 

2.023 

1.852 

1.670 

1.478 

1.282 

1.084 

0.893 
0.712 
0.548 
0.405 
0.285 
0.191 
0.121 
0.071 
0.039 
0 .020 
0.009 
o.oo4 
0.001 
0.001 
0 
0 


-2 .96 
-2.94 
-3-04 
-3 -24 
-3.53 
-3-90 
-4.27 
-4.79 
-5 .25 
-5 .67 
- 6 .02 
-6 .26 
-6.35 
-6 .28 
-6 .02 
-5*59 
- 5 .00 
-4.29 
-3.52 
-2 .74 
- 2.02 
-i.4o 
-0 .90 
-0 .54 

-0 .299 

-0 .150 
-0 .068 
-0 .027 
-0 .011 
-0 .004 
-0 .001 



I = 0.045 


0 

0 .0370 
0 .0745 
0 .1125 
0 .1509 
0 .1899 

0.2295 

0.2698 
0 .3107 
0 .3525 
0.3951 
0.4388 
0 .4837 
0.5297 
0 .5772 
0 .6262 
0 .6770 
0 .7297 
0 .7845 

0.8416 
0 .9011 
0.9633 
1 .0282 
1.0959 
1.1665 

1.2401 
1.3166 
1.3961 
1 .4786 
1 .5641 
1 .6526 


0 

0 .0805 
0 .1591 
0 .2356 
0.3098 
0.3817 
0 .4510 
0 .5174 
0.5805 
0 .64oo 
0.6956 
0 .7468 

0.7933 

0.8349 

0.8713 

0.9024 

0.9283 

0.9492 

0.9654 

0.9775 

0.9861 
0.9919 
0 .9956 
0 .9978 

0 .9990 

0 .9995 

0 .9998 
0 .9999 

1 
1 
1 


l 8u 

1 5 2 u 

u 0 3y* 

Uo Sy*2 

2.216 

-2.20 

2.136 

-2.15 

2.055 

-2.14 

1.973 

-2.19 

1 .888 

-2.27 

1.797 

-2.38 

1.701 

-2 .49 

1.597 

-2.67 

1.485 

-2.79 

1.365 

-2.91 

1.239 

-3 .00 

1.107 

-3-04 

0.972 

-3.01 

0.835 

-2.92 

0.700 

-2.75 

0.571 

-2.51 

0.451 

-2.21 

0.344 

-I.87 

O.252 

-1.51 

.0 .176 

-1.16 

0.116 

-0.84 

0.073 

-O.58 

0.043 

-O.367. 

0.023 

-0 .217 

0.012 

-0 .119 

0.005 

-0 .059 

0.002 

-0 .027 

0.001 

-0 .011 

0 

-0 .004 

0 

-0 .001 

0 

-0 .001 
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the velocity distribution with its first and 


SECOND DERIVATIVE FOR VARIOUS I 


cr = 0.4; | = 0.08 


a = 0.5; S = O.125 



1 chz 1 o^u 
Uo | Uo 3y* Uo 



1.835 

1.751 

1.669 

1.586 

1.502 

1.415 

1.326 

1.233 

1.136 

1.035 

0.931 

0.825 

0.717 

0.611 

0.508 

0.411 

O.322 

0.243 

0.176 

0.122 

0.080 

0.049 

0.029 

0.015 

0.008 

0.004 

0.001 

0.001 

0 

0 

0 


- 1.82 

-1.75 

-1.72 

-1.70 

-1.71 

-1.73 

-1.77 

- 1.80 

-1.84 

- 1.86 

-1.87 

-1.85 

-1.80 

-1.71 

-1.58 

-1.42 

-1.23 

- 1.02 

- 0.82 

- 0.62 

-0.444 

-0.299 

-0.186 

- 0.110 

-0.059 

-0.029 

- 0.012 

-0.005 

- 0.001 

0 

0 



1 du 1 ^u 

Uo 3F* U Q 2 


0.2314 
0.2921 
0.3542 
0.4176 
0.4826 
0.5492 
0 .'6176 
0.6880 
0.7605 
0.8354 
0.9129 
0.9933 
1.0768 
1.1638 
1.2545 
1.3493 

1.4482 
1.5516 
1.6596 
1.7724 
1.8900 
2 . 0126 
2.1401 
2.2726 
2.4101 
2.5526 
2.7001 


0.4073 

0.4783 

0.5454 

0.6084 

0.6671 

0.7211 

0.7704 

0.8145 

0.8534 

0.8871 

0.9154 

0.9387 

0.9572 

0.9713 

0.9816 

0.9889 

0.9937 

0.9966 

0.9983 

0.9992 

0.9997 

0.9999 

1.0000 

1 

1 

1 


1.612 

1.525 

1.439 

1.355 

1.271 

1.187 

1.101 

1.014 

0.926 

0.836 

0.745 

0.654 

0.564 

0.476 

0.393 

0.315 

0.244 

0.183 

O.131 

0.090 

0.058 

0.036 

0.020 


-1.60 
-1.52 
-1.47 
-1.43 
-1.40 
-1.38 
-1.37 
-1.37 
-1.36 
-1.34 
-1.31 
-1.27 
- 1.21 
-1.13 
-1.03 
-0.91 
-0.78 
-0.64 
-0.50 
-0.376 
— 0.266 
-0.177 
- 0.110 
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TABLE I.— Continued 


TEE "TELOCITY DISTRIBUTION WITH ITS 


SECOND DERIVATIVE FOR VARIOUS 



VO 

. 

0 

II 

to 

w * 

II 

O 

. 

H 

CD 

a 

~ V 0 J 

y* - ° 
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u 

Do 

1 5u 
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1 ' d 2 u 

Do dy*2 

y* = " T ° y 
J V 

0 

0 

1.467 

-1.45 

0 

0.06I7 

0.0919 

1.376 

-1.37 

0.0724 

0.1307 

0.1798 

1.287 

-1.30 

0.1464 

0.1982 

0.2637 

■ 1.201 

-1.25 

0.2225 

0.2672 

0.3437 

1.117 

-1.20 


0.3379 

0.4197 

1.033 

-1.16 

0.3805 

0.4102 

0.4949 

0.950 

-1.13 

0.4627 

0.4844 

0.5593 

0.868 

-1.101 

0.5471 

0.5606 

0.6218 

0.785 

-1.07 

0.6341 

0.6388 

0.6800 

0.703 

-1.04 

0.7237 

0.7194 

0.7333 

0.621 

-0.99 

O.8161 

0.8026 

0.7815 

0.540 

-0.94 

0.9117 

0.8884 

0.8245 

0.462 

— 0. 88 

1.0106 

0.9773 

0.8621 

0.386 

-0.81 

.1.1132 

1.0695 

0.8944 

0.316 

-0.73 

1.2206 

1.1653 

0.9215 

0.251 

—0. 63 

1.3316 

1.2649 

0.9435 

0.193 

-^.532 

1.4473 

1.3689 

0.9608 

0.143 

-0.413 

1.5680 

1.4774 

0.9735 

0.101 

-0.334 

1.6942 

1.5907 

0.9835 

0.069 

-0.247 

1.8262 

1.7093 

0.9901 

0.044 

-0.172 

1.9643 

1.8332 

0.9944 

0.027 

— 0.113 

2.1088 

1.9628 

0.9970 

0.015 

— 0. 069 

2.2599 

2.0981 

O.9985 

0.008 

-0.039 

2.4177 

I 2.2392 

0.9993 

0.004 

-0.020 

2.5823 

| 2.3863 

0.9997 

0.002 

-0.010 

2.7539 

! 2.5393 

0.9999 

0.001 

-0.004 

2.9324 

1 2.6983 

0.9999 

0 

-0.002 

3. 1179 


1 

0 

-0.001 

3.3104 

3.0343 

1 

0 

0 

3.5099 

3.2113 

1 



0 

0 

3.7164 


NACA TM No. 1205 


FIRST AND 
i 


11 

0 
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-'■a 

L/TT 

= 0.245 


u 

1 8u 

1 8 2 u 

Do 

Do by * 

Do Sy* 2 

0 

I.366 

-1.35 

0.0955 

I.271 

-1.27 

0.1864 

1.180 

-1.19 

O.2727 

I.O92 

-1.12 

0.3545 

1.007 

-1.07 

0.4317 

O.923 

-1.02 

0.5042 

0.842 

-0.97 

0.5719 

0.762 

-0.93 

O.6347 

O.683 

-0.88 * 

0.6924 

O.606 

-0.84 

0.7450 

0.531 

-0.79 

0.7922 

0.458 

-0.74 

0.8340 

O.388 

— 0. 68 

O.8704 

0.322 

-O.61 

O.9013 

0.261 

-0.54 

0.9270 

0.205 

-0.462 

0.9479 

O.I56 

-0.384 

0.9642 

0.115 

-0.307 

0.9764 

0.081 

-O.235 

O.9852 

0.054 

-O.171 

O.9907 

0.034 

-0.117 

0.9951 

0.020 

-0.076 

0.9974 

0.011 

-0.046 

0.9987 

0.006 

-0.026 

0.9994 

0.003 

-O.013 

0.9998 

0.001 

-0.006 

0.9999 

0.001 

-0.003 

1 

0 

-0.001 

1 

0 

-0.001 

1 

0 

0 

1 

0 

0 
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0 

O.OT97 

0.l6l6 

0.2457 

O.3321 

0.4211 

0.5126 

0.6070 

0.7043 

0.8048 

0.9088 

I.OI65 

1.1283 

-1 ol.l.i, 


1.3652 

1.4912 

1.6228 

1.7603 

1.9041 

2.0547 

2.2123 

2.3773 

2.5498 

2.7301 

2.9183 

3.1143 

3.3184 

3.5304 

3.7504 

3.9784 

4.2144 


0 

0.0990 

0.1929 

0.2815 

O.365I 

0.4434 

0.5166 

0.5846 

0.6472 

0.7044 

0.7562 

0.8024 

0.8431 

0.8782 

O.9079 

0.9324 

O.9521 

0.9673 

0.9786 

0.9867 

0.9922 

0.9957 

0.9977 

0.9989 

0.9995 

0.9998 

0.9999 

1 

l 

1 

1 


1.292 

1.194 

1.100 

1.010 

0.923 

0.840 

0.760 

0.682 

0.606 

0.533 

0.463 

0.396 

0.333 

0.274 

0.219 

0.171 

0.129 

0.094 

0.065 

0.043 

0.027 

0.016 

0.009 

0.005 

0.002 

0.001 

0 

0 

0 

0 


-1.03 

-0.97 

-0.91 

-0.85 

-0.80 

-0.75 

-0.70 

-0.65 

-0.60 

-0.54 

-0.479 

-0.416 

-0.352 

- 0.288 

-0.227 

-0.171 

-0.123 

-0.083 

-0.053 

-0.032 

-0.018 

-0.009 

-0.004 

- 0.002 

- 0.001 

0 

0 


0 

0.0865 

0.1756 

0.2673 

0.3618 

0.4592 

0.5597 

0.6635 

0.7709 

0.8820 

0.9972 

1.1168 

1.2411 

1.3706 

1.5055 

1.6464 

1.7937 

1.9478 

2.1093 

2.2784 

2.4555 

2.6410 

2.8350 

3.0378 

3,2494 

3.4700 

3.6995 

3.9380 

4.1855 

4.4420 

4.7075 


0.8517 

0.8856 

0.9141 

0.9374 

0.9559 

O.9701 

0.9807 

0.9881 

0.9930 

0.9962 

0.9980 


0.347 

0.289 

0.236 

0.187 

0.145 

0.108 

0.078 

0.054. 

0.035 

0.022 

0.013 

0.007 

0.004 

0.002 



-O.96 

-0.89 

-O.83 

-0.76 

-0.71 

-O.65 

-0.60 

-0.55 

-0.495 

-0.441 

-O.386 

-0.331 

-O.276 

- 0.222 

-0.173 

-O.129 

-O.O9I 

-O.O6I 

-O.O38 

-0.023 

- 0.012 

-O.OO6 

-O.OO3 

- 0.001 

0 

0 

0 

0 
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~ Y 0 y 

V 


0 

0.0929 

0.1888 

0.2877 

O.3898 

0.4953 

0.6044 

0.7174 

O.8345 

O.9560 


T = 1.0; 

i = 0.5 


u 

1 8 u 

1 8 ^ . 

% 

Oo 8 y* 

3y*2 


0 

0.1060 

0.2055 

0.2986 

O.3854 

0.4660 

0.5405 

0.6088 

0.6710 

0.7272 


1.194 

i.089 

0.989 

0.895 

0.806 

0.722 

0.643 

O.568 

0.497 

0.430 

0.367 

0.308 

0.254 

0.206 

0.162 

0.124 



0.001 

0 

0 

0 

0 

0 


-1.18 
-1.08 
-0.99 
-0.91 
-0.83 
-0.76 
-0.70 
— 0 . 64 
-0.58 
-0.52 
-0.472 

- 0.420 

-O.369 

-0.318 

- 0.268 

- 0.221 

- 0.176 

-O.I35- 

- 0.099 

-O.O69 

- 0.046 

- 0.028 

-0.016 

-0.009 

- 0.004 

- 0.002 

- 0.001 

0 

0 

0 

0 



0 

0.1047 

0.2130 

0.3253 

0.4417 

0.5625 

0.6880 

0.8184 

0.9541 

1.0955 

1.2430 

1.3969 

1.5578 

1.7261 

1.9025 

2.0874 

2.2814 

2.4851 
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Figure 4.- Sketch of coordinates. 
Vr 



Figure 6. - Sketch of coordinates. 
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Figure 8.- 


Velocity distribution 










64 


NACA TM No. 1205 



Figure 9. - The first derivative of the velocity distribution 
6*/U 0 5u/dy against y/8*. 











Figure 15.- The universal law for the friction drag of the plate with 
homogeneous suction F(|j) = c f /2 c against g = c 2 U l A> 












